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14- Abstroct 

Experxmental values of acoustic wave propagation constant and 
characteristic impedance in fibrous materials, and normal' absorption for two 
plastic foams, have been compared to theoretical predictions obtained with Biot’s 
theory. The best agreement has been observed- for. fibrous materials between 
Biot’s theory and Delany and Bazley experiments for a nearly zero mass coupling 
parameter. For foams, the ^/4 structure resonance effect on absorption has 
been calculated by usipg four-pole modelling of the medium. A significant mass 
coupling parameter. is then necessary for obtaining agreement between the 
behavior of the measured absorption coefficients and the theoretical predictions. 
It is shown how the formalism used for predicting foams absorption coefficients 
may be used for* s-tudying the acoustic behavior of multi-layered media. 
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BIOT THEORY AND ACOUSTICAL PROPERTIES OF HIGH POROSITY 
FIBROUS MATERIALS AND PLASTIC FOAMS 

J. Allard and A. Aknine 

Acoustical Laboratory, Faculte des Sciences du Mans 
1. Introduction /221 * 

On an experimental level, the absorption and propagation of 
sound waves in fibrous materials are well described by the 
empirical laws of Delany and Bazley [1, 2], However, these laws 
do not account for the effect of vibration of the fibrous 
structure when it is significant, and are only valid in a 
limited frequency range. 

The acoustical properties of plastic foams are much /222 

less well known; the only systematic study comparable to that of 
Delany and Bazley, done by Cummings [3], has bearing on only one 
type of material. On a theoretical level, there are numerous 
studies. Attenborough did a very complete review of work done 
up to 1981 [4]. 

We refer in particular to studies done by Sides, 
Attenborough, and Mullholland [5] and by Moore and Lyon [6], 
using the theories of Zwikker and Kosten [7] and of Biot C8, 9] 
with some modifications. 

Sides [5] uses Biot's theory [8] without the mass coupling 
parameter and comes up with a good approximation of Delany' s and 
Bazley 's laws Cl]. Moore [6] brings out analogies between 
Biot's theory [8] and that of Zwikker and Kosten [7] and uses a 
description of the porous materials which is quite close to 
these two theories; however, no reference is made to 
experimental results. 


♦Numbers in the margin indicate pagination in the foreign text. 
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We have developed a general formalism based on Biot's 
theory. This formalism allows for transmission of vibrations by 
a material used for thermal insulation, consisting of a glass 
wool layer between two layers of plaster. It also allows for 
sonic absorption by a panel installed on a rigid surface. 

We have analyzed experimental results on fibrous materials 
and two plastic foams, and it was possible to reveal the effect 
of the mass coupling parameter and quarter-wave resonances of 
the structure. 


2. Description of the propagation of sound in high-porosity 
materials according to Biot’s theory 

2.1. Stress-deformation equations 

The material is assumed to be of elastic structure. It is 
also assumed to be homogeneous and isotropic on the wavelength 
scale. The quantities below are meaningful only in this 
macroscopic context. 


Let o\ j and e^^ be the components of the stress and 
deformation tensors for the structure; e and e the expansions of 
the structure and of the air in the material. 


Stresses are the forces per unit of surface of material 
acting on the structure. 

Making p^ the pressure in the fluid, it is possible also 
to define a normal force per unit of surface of material acting 
on the gas; 



( 1 ) 


h is the porosity of the material 



From Biot's energy considerations [8], the following 
equations between normal stresses and deformations were obtained 

o xx = 2 N e xx A- A e 4- Q c 
<j vv =“lN€ vv +Ae-\- Q c 

( 2 ) 

°z z =2Ne :: +AeA-Qe ' ' 

$— Qe+ e R 

where N is the shearing module of the structure • Biot and 
Willis CIO] proposed two experiments to measure two intermediate 
quantities for evaluation of A, Q, and R. 

The first quantity is the compressibility yc of a sample of 
fibrous or porous material enveloped by a supple membrane 
pierced by a hole so that the pressure in the sample remains 
constant and equal to the exterior pressure during action on the 
envelope . 

The second quantity is the compressibility 6 of an 
unenveloped sample during variation of exterior pressure, i.e. 
the pressure of the material making up the structure, e.g. that 
of plain glass for glass wool. The experiments proposed by Biot 
are very symbolic in a dynamic regime, but the representation 
of x and b remains as stated above. 


For high-porosity materials like most glass or stone wools, 
or the foams studied for which bubble walls are in large part 
destroyed, <5 is negligeable before yt and compressibility 1/y of 
the air, and it can be stated that: 


5 = o. 


( 3 ) 


On the other hand, it can be assumed, as was shown for 
fibrous materials [4], that the Poisson coefficient v of the 
structure is zero for all the high-porosity materials studied. 



It is therefore possible to write: 


1/H = 2N/3. (4) 

The dynamic shearing module N can be measured at atmospheric 
pressure, because the effect of the air in the porous material 
is only a slight increase in the weight of the structure. 

The general equations expressing A, Q, and R as a function 
of h, N, ~K, and & [10] are simplified if h is made equal to 1 
and if equations (3) and (4) are verified. The following is 
obtained: 


A — 0, Q - 0, R = l/y * 

(5) 

According to [7], ^ can be written: 
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/ 0,8 y,(0,86/^\ 

* ’• ° \ 0,86 nfjh (0,86// 1/7/ ' 


In this expression, is the quantity: 


fi = y$O)kQ 0 fh o . 

(6) 


a 


Q o 

k 


is the specific resistance of the material to the passage 
of air. 

is the pulsation of the sound wave, 
is the volumic mass of the air. 

is the shape factor from Zwikker and Kosten [7]. 


This factor, for a material with parallel cylindrical pores 
of the same diameter, is the inverse of the square of the cosine 
of the pore generator angle with normal on the surface. It is 
associated with the average inclination of the airflow with 
respect to normal on the surface. 



2.2. Equations for movement and quantities characterizing the 
propagation of sound waves 

Let : 

be the structure mass per unit of volume of material. 

£>2 be the air mass per unit of volume of material 

^2 = h ^0' 

£ a be the mass coupling parameter introduced by Biot [8] based 
on energy considerations. 

£?l 1 = (?l + 6a 

On = Qi + 0a 
(?12 = — • 

The paired equations for the expansion waves in Biot's theory 
[8] are expressed: 

02 0 

V 2 (Pe+ Qe)=j^(eu +e\2 e) + b-(e-e) 

V J (0 e + R e) = (g n e + 622*) ~ bjj{e - c) . ^ 8 ^ 

8 

The quantity b — {e — e) is a damping parameter related to the 
viscosity of the air. 


By analogy with the porous materials with cylindrical 
pores, it is possible to write [5, 7]: 


I^T J\ (p VJ ) 

g ^ 7 Join ]fj) 

4 1 2y » k. 

v YJ Jo (/« VJ) 


Moore [6] remarked that Q can be identified with the quantity 

cl 

(k - 1)6 q which appears in Zwikker and Kosten [7] (page 53): 


n 


(k - 1 )£ 0 . 


( 9 ) 


From equation 8, the two propagation constants 
obtained from the expansion waves by writing: 


and y b are 


e — a l exp j (yx - co t) 

£ = a 2 t\ pj (yx — cot) . 


For each of these values of y, it is possible to obtain with one 
or the other of equations (3) the ratios and <£> b between 
the expansion amplitudes in the fluid and the solid: 

f- a 2 / a i* 

Finally, for each of the progressive waves associated with y a 
and y^, two characteristic impedances and are 
defined, which are the ratio of the forces per unit of surface 
for the structure and the gas to the velocities of the structure 
and of the gas. 

In the case of a planar longitudinal wave propagating in 
direction Ox, the following is obtained from equations (2) and 
(5): 


Z] = 2NyJco 
Z^lN^/co 
Z\ = R yjco 
Z^R^/oj. 


( 10 ) 


3. Empirical laws of Delany and Bazley and Biot's theory - 
evaluation of the mass coupling parameter for fibrous 
materials 

3.1. Qualitative aspect 

The propagation of sound in fibrous materials for 
frequencies f varying between 0.01 a and C can be associated 
with a wave in the gas of the material whose characteristic 
impedance and propagation constant depend only on f/a, according 
to the empirical laws of Delany and Bazley Cl], We have 
verified the compatibility of these empirical laws with Biot's 




theory [8] for a common type of glass wool with porosity h near 
1# Poisson coefficient v near 0, volumic mass q , and specific 
resistance to the passage of air 0 with values: 

a i = 30 kg/m 3 
cr = 20,000 kg m -3 s" 1 . 

The empirical laws of Delany and Bazley [1] for this material 

are thus verified for f varying between 200 Hz and 20,000 Hz. 

For N varying between 10 and 10 N/m and for frequencies 

above 200 Hz, one of the two propagation constants from /224 

Biot's theiry, written as is such that J(Pgf is between 4 

and 20. The associated wave thus propagates essentially in the 

gas. The characteristic impedance associated with the 

propagation of this wave in the gas is written as zf[. The 

other wave has a propagation constant of y and a 

s ' s 

characteristic impedance of Z^, which are approximately 
those of the structure in a vacuum. 



( I . -3 

<p s I varies between 1 and 10 . The structure thus pushes 

against air as it moves, but since it is much heavier than air, 

this does not noticeably change its movement. 

This phenomenon of partial decoupling of the structure was 
brought to light by Zwikker and Kosten [7] and exists only for f 
verifying: 

C7/2TT f 1, (12) 

or, for the material considered, f >100 Hz. 

For greater values of cr/ f, the preceding description is no 
longer valid. 


<7 



The contribution of this wave is slight in experiments to 
measure impedances, because it is difficult for the air 
surrounding the structure to move it, except in the cases of 
quarter-wave resonances of the structure, which will be examined 
later. This means that Delany and Bazley [1] did not reveal two 
waves of different velocities or phenomena linked to vibrations 
of the structure. 


3.2. 


Comparison of the values of and 2^ obtained 

using Biot's theory with those obtained using empirical 

laws and determination of the mass coupling parameter 


The quantities y^ and 
rigidity of the structure. 


vary very little with the 
as shown in Table 1 . 


3 2 

We have thus arbitrarily written N = 10 N/m without 
this noticeably influencing later results, and, for the fibrous 
material considered, we have compared the predictions from 
Biot's theory [8] with the empirical laws of Delany and Bazley 
[1]. The following quantities are defined with Delany* s and 
Bazley 's notations: 


CL = Im (y ). R = Re(z|), 

/$ = Re(v ), X = Im(z|). 

These quantities depend on form factor k. 

In the validity zone of the empirical laws, values of f /o 
were sought for which these quantities vary significantly as a 
function of k. Only the quantities and R vary noticeably, and 
this for high values of f /<j . Putting aside f/cr = 1, which is at 
the extreme end of the validity zone, in Table II we have shown 
the variations of p and R as a function of k for f/<j = 0.1 and 
f/o = 0.316. The last column represents the values obtained 
using empirical laws. The speed of sound in free air is written 
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Table I. 


Variations of V g and Z§ as a function of N 
' and k = 1.3. 


for f = 800 Hz 


Shearing 

Module 

N 

M/m^ 


Propagation 

Constant 



Characteristic 

Impedance 

Z§J 

Nm - 3 S 


10 s 0,30- 10 2 +j0,17- 10 2 

10 3 0,29 - 1 0 2 + j 0. 1 6 • 10 J 


0,64- 10 3 +j 0,29 10 3 
0,63- l0 3 +j0,29- 10 3 


Table II. 

Variation of the quantities jS and R as a function of k for 

t/a = 0.1 and 0.316. 


k 

Co/? _ | 
to 

ior//(7= 

R i 

ft i Co 

0,316kg-'m 3 

Cp/? | 

CO 

$or f/o — 

R { 

Qo Q 

0,1 kg~ 1 m 3 

1 

0,24 

0,08 

0,45 

0,16 

1,1 

0,29 

0,15 

0,49 

0,22 

1,2 

0,34 

0,21 

0,53 

0,27 

1,3 

0,39 

0,26 

0,58 

0,33 

1,4 

0,43 

0,30 

0,62 

0,39 

1,5 

0,48 

0,35 

0,66 

0,44 

Delany et 
Bazley 

0,20 

0,12 

0,43 

0.29 


Table III. 

Values of d, , R, and X obtained with the laws of Delany and 
Bazley (D.B.) and Biot's theory for a 0 mass coupling parameter. 


/ 

C 0 a 

Co/? 1 

w 

R 

. * 
Qo Q) 

Ref. 

o 

kg“' m 3 

CO 

Oo Q 

0,0316 

1,34 

1,4 

0,96 

u 

0,68 

0,56 

0,96 

0,92 

D.B. 

Biol 

0,1 

0,68 

0,67 

0,43 

0,45 

0,29 

0,16 

0,41 

0,40 

D.B 

Biot 

0,316 

0,35 

0,30 

0,20 

0,24 

0,12 

0,08 

0,18 

0,14 

D.B- 

Biot 


For f/cr = 0.316 and 0.1, the best agreement with 
experimental results is obtained for values of k between 1 and 

1.2, i.e., according to equation (9), for values of g between 

3 a 

0 and 0.25 kg/m . In fact, a mass coupling parameter other 
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than 0 does not noticeably improve agreement with experi- 
mental results, and for mass coupling parameters greater than 
0.25 kg/m , agreement worsens. The mass coupling parameter is 
thus negligeable for fibrous materials. 

For £> a = 0, Table III shows the values of cL, R, and X 
obtained using Biot's theory [8] and Delany's and Bazley's laws 
[1] within their validity range. Agreement is generally 
excellent . 


4. Modelling of longitudinally vibrating porous materials 
4.1. Chain matrix and equivalent four-pole 

A thickness L of material will be represented by the four- 
pole in Figure 1. 

The chain matrix of the four-pole makes it possible to 
determine forces p^ and p 2 per unit of surface on the 
structure and the gas, as well as velocities v^ and of 
the structure and the gas on the front face using the same 
quantities defined on the back face. 



We write: 

4,-Zt Z\<p*-Z\Z\<Pi 

A 1 =Z\Z\<p i ~Z\<p^Z 3 i . 
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The elements of the chain matrix are expressed as follows: 

Ail = (Zf Z\ <p b cos y t L — Z\ Z\ (p a cos y b L)/A , ! 

(cos y a L - cos y b £)/Z , 

Aji = j (- Z\(p h sin y a L + Z\ <p a sin >- b L)/A , 

Aii = j (- <p a Z\ sin y a L+<p h Z\ <p z sin y b L)/A , 

A12 = (Z? Z f) (cos y a L- cos y b L)/A 2 I 

Aa = (Z| <p 3 Z \ cos >’ a L- Z: f b Z, cos y b L)/J 2 
A2 = j (— Z ^ 1 sin'/ a L+ Z] siny b L)/A 2 
6 « = j (~ (ffaZ^sin/a L+ <hZ]siny b L)/A2 
An = j (- Z^ b sin y a L + Z\tp a sin y b L)/(^> b - <p a ) 

V 3 ~ J (- Z a 2 <p a p b sin y a L + Z£? a sin y b L)/(«> b - ? a ) 

-An = (<o b cos y a L- <p. d cos y b L)/(<p b - ^ a ) | 

h ^ ^ (cos y a A - cos y b £)/(« 5 b - </? a ) 

|V. j 

jV® ffl T ( Z ' Sitl >,b L “ Z “ Sin Za Z-) 

j (Z^ p b sin y b L- z\(p. i sin y a L)/(<p a - ^ b ) 

§* * (cos >' a I - cos y b L)/(<p a - <p h ) 

§* ~ (^ cos y a L - <p b cos y b L)/(<p a - p b ) . 


HIM 

iL jp? 


® ® 


Figure 1. 

a) Sample with thickness L; b) Equivalent four-pole 



4.2. Chain matrix associated with a group of several 
superimposed layers of porous materials 

The chain matrix allows easy description of the 
superimposed layers of porous materials. 

Assume two layers A and B. At the interface of these two 
materials, continuity equations (14) must be verified: 




P 2A^ h A P 2B/ h B 
P 1A + P 2A = P 1B + P 2B 


h A (v 2a 


" V 1A> = 


h B (v 2B 


Continuity of the structure's 
velocity 

Pressure continuity (14) 
Continuity of the total 
force per unit of surface 
v 1 b) Continuity of air flux 


We write [b A ] and [& B ] for chain matrices A and B and [S ] 
for the interface matrix which, according to equations (14), 
expressed: 


"I 1 -A a //i b 

_ 0 /i A //i B 

ldl3 0 0 

.0 0 


0 0 ■ 

0 0 

1 0 

l-h b /h A V/»a. 


is 


It is the unit matrix if environments A and B have the same 
porosity. 

The product [£ A ] C U CS>g] is the chain matrix of 
group AB. 


5. Effect of the structure's quarter-wave resonance on the 

absorption of a plastic foam panel placed on a rigid support 

The presence of the structure's quarter-wave resonance can 
considerably modify absorption with respect to the predictions 
of Delany and Bazley [l]y the wave associated with the 
structure's vibrations can thus be of high amplitude. 

For many plastic foam panels, this resonance, calculated 
from equation (11), is between 500 Hz and 1500 Hz. 


Absorption will be obtained from impedance Z on the /226 

surface of the panel calculated by writing, in equation (13): 


= v„ = 


/y 


0. 


(15) 



For materials with porosity near 1, it is also possible to write: 


0 0 0 

V 2 = v ' p 2 = p ' P 1 = °* (16) 

p and v are the pressure and the acoustical velocity at the 
surface of the sample. 


Using equations (13), (15), and (16), Z is deduced from the 
elements of the chain matrix and normal incidence absorption. 



Figure 2. 

Variations of the absorption coefficient of foam 1 as a function 
of frequency f. 

Key : A - Absorption coefficient; B - Frequency. 

Curve a: Experimentation, 

Curve b: k = 1, zero mass coupling coefficient (Biot theory). 
Curve c: k = 3, mass coupling coefficient equal to twice the 
volumic mass of the air (Biot theory). 



Figure 3. 

Variations of the absorption coefficient of foam 2 as a function 
of frequency f. 

Key : A - Absorption coefficient; B - Frequency. 

Curve a: Experimentation, 

Curve b: k = 1, zero mass coupling coefficient (Biot theory). 
Curve c: k = 3, mass coupling coefficient equal to twice the 
volumic mass of the air (Biot theory). 




Absorption was calculated as a function of frequency for 
two plastic foams manufactured by Roth Freres (Roth Freres S.A., 
BP 13, 67023 Strasbourg Cedex) with porosity near 1 and Poisson 
coefficient assumed to be near 0. 

The other characteristics of these two materials are given 
in Table IV. 


Table IV. 

Characteristics of foams 1 and 2. 


Foam Thickness 


L 

m 


Volumic 
mass of 
the 

structure 


m^ 


Specific 
resistance 
to the 
passage 
of air 

kg m" J s~ x 


Dynamic 

shearing 

module 


N 

N/m 2 



1 4. 8*10- 2 9.7 7.7-10 3 1.4-10 5 

2 3.7-10- 2 18 16-10 3 1.4-10 5 


The normal absorption calculated according to the preceding 
method for k = 1 and k = 3 and measured experimentally is shown 
as a function of frequency in Figure 2 and Figure 3. There is a 
certain dispersion in the experimental curves, because the 
measurements were made in the Kundt tube, and this method is not 
well adapted for supple materials [113 if the vibrations of the 
structure participate in absorption. 

The quarter wave resonances for the two samples are 
localized respectively at 890 Hz and 850 Hz. 

In the area of these resonances, the absorption 
coefficients measured are very different from those obtained 
with Delany's and Bazley's empirical laws [1]. 

With the calculation method used, the agreement between the 
experimental and calculated absorption coefficients is not good 
for k = 1, but it is much better for k = 3, i.e. when a 



significant mass coupling parameter Q. a is used, on the order 
of twice the volumic mass of the air. 


6. Conclusion 

The Biot theory makes it possible to correctly predict the 
acoustical behavior of fibrous materials with a negligeable mass 
coupling parameter. Regarding plastic foams, no experimental 
study has been done on a large number of materials and under 
good conditions, i.e. in a free field on large-sized samples 
such that the structure can vibrate freely under conditions 
similar to those under which the material is used. 

However, for the two samples that were studied, the / 277 

Biot theory made it possible to predict absorption behavior as a 
function of frequency when a mass coupling parameter on the 
order of twice the volumic mass of the air was used. 

The trajectories of air molecules around the walls of the 
bubbles in foams can be much more perturbed than in fibrous 
materials, which explains the importance of form and thus of 
mass coupling in foams. 

Finally, note the remarkable variety of areas where the 
Biot theory, initially conceived to describe the behavior of 
rocks saturated with liquid, can be used. In addition to its 
applications in the area of sound waves, this theory was found 
to be useful in describing the propagation of ultrasonic waves, 
for example in glass spheres saturated with water [12]. 

(Received April 27, 1984) 
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